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ABSTRACT 

We present a physically motivated model to estimate the molecular hydrogen (H2) 
content of high-redshift (z w 5.7, 6.6) Lyman Alpha Emitters (LAEs) extracted from 
a suite of cosmological simulations. We find that the H2 mass fraction, /h 2 , depends 
on three main LAE physical properties: (a) star formation rate, (b) dust mass, and 
(c) cold neutral gas mass. At z 5.7, the value of / H2 peaks and ranges between 
0.5 — 0.9 for intermediate mass LAEs with stellar mass M* ss 10 9 ~ M©, decreasing 
for both smaller and larger galaxies. However, the largest value of the H2 mass is 
found in the most luminous LAEs. These trends also hold at z « 6.6, although, due 
to a lower dust content, fn 2 (z = 6.6) ~ 0.5/h 2 (z = 5.7) when averaged over all LAEs; 
they arise due to the interplay between the H2 formation/shielding controlled by dust 
and the intensity of the ultraviolet (UV) Lyman- Werner photo-dissociating radiation 
produced by stars. We then predict the carbon monoxide (CO) luminosities for such 
LAEs and check that they are consistent with the upper limits found by Wagg et 
al. (2009) for two z > 6 LAEs. At z 5.7,6.6, the lowest CO rotational transition 
observable for both samples with the actual capabilities of Atacama Large Millimeter 
Array (ALMA) is the CO(6-5). We find that atz« 5.7, about 1-2% of LAEs, i.e., 
those with an observed Lyman Alpha luminosity larger than 10 43 2 ergs _1 would be 
detectable with an integration time of 5-10 hours (S/N = 5); at z f=s 6.6 none of 
the LAEs would be detectable in CO, even with an ALMA integration time of 10 
hours. We also build the CO 'flux function', i.e., the number density of LAEs as a 
function of the line-integrated CO flux, S'coj an d show that it peaks at Sqo = 0.1 
mJy at z = 5.7, progressively shifting to lower values at higher redshifts. We end by 
discussing the model uncertainties. 
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radio lines: ISM, cosmology: theory 



1 INTRODUCTION 

The past few years have seen a rapid increase in the amount 
of data available on high-redshift galaxies. This has been 
made possible by a combination of state of the art instru- 
ments such as the Hubble Space Telescope, the Subaru and 
Keck telescopes, and sophisticated selection techniques. Of 
the latter, one of the most successful approaches has been 
the use of the narrow-band technique (e.g. Malhotra et al. 
20051 JStumasaku et al.] [20061 |Kashikawa et al.||2006|lHu 



et al.| 20101 that is based on looking for the Lyman al- 
pha (Lya) emission at 1216 A in the galaxy rest frame. 
Hundreds of such galaxies, called LAEs have now been con- 
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7 ( |Iye et al.|2006[ ).~Due to their large number statistics 
and unambiguous spectral signatures, LAEs are arguably 
the best probes of reionization and high-redshift galaxy evo- 
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However, using LAEs as probes of high-redshift galaxy 
populations, understanding them to study reionization, and 
calculating their contribution to reionization requires an un- 
derstanding of their star formation rates (SFR). Translat- 
ing the observed UV luminosity (1375 A in the galaxy rest 
frame) into an intrinsic SFR is rendered hard by the fact that 
the continuum photons produced in a galaxy are attenuated 
by the dust in the galactic interstellar medium (ISM) before 
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they reach the observer. Inferring the intrinsic SFR using 
the observed Lya luminosity is even more complicated since 
Lya photons are absorbed both by the ISM dust, as well 
as the neutral hydrogen (H i ) in the intergalactic medium 
(IGM) along the line of sight between the emitter and the 
observer. Indeed, a number of studies point to LAEs being 
dust enriched, even at redshifts z ~ 6: using theoretical mod- 
els, Dayal et al. ( 2010 1 have shown that at z ~ 5.7, the color 
excess of LAEs, E(B — V) ~ 0.15 while observationally, the 
color excess values range between E(B — V) ~ 0.025 — 0.32 
at z ss 4 - 5.7 dLai et al.||2007| |Pirzkal et aT||2007| [Flnkel- 
stein et al.|2009 1. Further, the observed Lya luminosity de- 
pends both on the reionization state of the IGM, as well as 
on the IGM peculiar velocities along the line of sight; in- 
flows/outflows into/from a galaxy can blueshift/redshift the 
Lya line, thereby decreasing/increasing the IGM Lya trans- 



mission, T a (|Verhamme et al.|2006] |Dayal fc Ferrara|2011b 



Dijkstra et al. 20111. However, the extent to which pecu- 



liar velocities influence T a is debatable since these calcula- 
tions have mostly been performed under idealized situations. 



For example, Verhamme et al. ( 2006 \ have used spherically 



symmetric outflows of H i to show an enhancement in T a 



however, many studies, e.g. Fangano et al. (2007) and refer- 



ences therein, have shown that Kelvin-Helmholtz instabili- 
ties would result in breaking-up such symmetric outflows. 

In this sense, H2 is a far better indicator of the SFR 
since stars form in dense, cold molecular clouds (MCs); the- 
oretical and observational constraints on the latter are then 
of utmost importance to shed light on the molecular content, 
and therefore the intrinsic SFR of these high-redshift galax- 
ies. The H2 content of galaxies is generally studied through 



observations of CO rotational emission lines (Solomon & 



Vanden Bout 2005 Omont 20071 that have been detected 



in more than a hundred high-redshift sources, even though 
searches for molecular gas at redshifts z > 4 have so far 
focused mainly on quasars and on the most massive, far- 



infrared-luminous and submillimeter galaxies (see Riechers 



2011 1. As of now, only scant effort has been devoted to ob- 



serving the molecular content of high-redshift (z ^ 6) LAEs. 
In one of the only observational works available, Wagg et al. 
( 2009 1 have searched for low-J rotational CO emission lines 



in two LAEs at z > 6.5; the non-detection of any CO emis- 
sion from these galaxies can then be used to put constraints 
on the amount of molecular gas in these sources. 

In this work, our aim is to present a self-consistent and 
physically motivated model to calculate the H2 fraction and 
mass, and relate it to the physical properties of the galaxies 
identified as LAEs at z « 5.7, 6.6. To do so, we start by using 
a previously developed LAE model (see |Dayal et al.||2008| 
2009 2010; Da yal fc Ferrara|2011a I, where the authors com- 
bined state of the art cosmological SPH simulations with a 
Lya production/transmission model to successfully repro- 
duce a large number of observational data sets. We couple 
this with a semi-analytic model that describes the structure 
of the MCs, to calculate the molecular hydrogen content 
taking into account its formation on dust grains, its destruc- 
tion by UV photons, and the shielding by H 1 in the ISM, 
as proposed by|Krumholz et al.| ( |2008| |2009| ) and |McKee~fc| 
Krumholz (2010). Once the molecular fraction is calculated 
for all of the LAEs at z ~ 5.7,6.6, we examine its correla- 
tions with the physical properties of the emitters, including 
the SFR, the total ISM dust mass and the amount of cold 



H 1 gas in the ISM. Translating the total H2 mass into a 
CO luminosity, we compare the predictions of our model to 



the observations of Wagg et al. (20091. Finally, we estimate 
the time required to observe such CO emission with ALMA, 
one of whose main goals is to observe molecular gas in high-z 
sources. 



2 THE COSMOLOGICAL SIMULATIONS 

In this Section, we briefly describe the simulation used in 
this work and interested readers are referred to ITornatorel 



et al. (2010) for a complete description. The simulation has 



been carried out using the TreePM-SPH code GADGET-2 



(Springcl 2005) with the implementation of chemodynam- 
ics as described in Tornatore et al. (20071. The adopted 



cosmological model corresponds to a ACDM Universe with 
Qm = 0.26, n A = 0.74, Q b = 0.0413, n s = 0.95, H = 
73kms _1 Mpc -1 and ag = 0.8, thus consistent with the 5- 



year analysis of the WMAP data (Komatsu et al. 20091 



The simulation has a periodic box size of 75h~ comoving 
Mpc (cMpc) and contains 512 3 dark matter (DM) particles, 
and initially the same number of baryonic particles. The 
run assumes a uniform redshift-dependent UV background 
produced by quasars and galaxies, as given by |Haardt "fc| 
|Madau| ( |1996[) and includes me tallicity-dependent radiative 
cooling (Sutherland & Dopita 1993|. The code has an ef- 
fective model to describe star formation from a multi-phase 
ISM and a prescription for galactic winds triggered by su- 



pernova (SN) explosions (Springel & Hernquist 20031; the 
initial mass function (IMF) is taken to be Salpeter between 
1 — 100 Mq. Metals and energy are released by stars of dif- 
ferent masses by properly accounting for mass-dependent 



lifetimes as proposed by Padovani & Matteucci (19931. The 



code uses the metallicity-dependent yields from 



Woosley &q 



Weaver ( |1995| ; the yields for SNIa and asymptotic giant 



branch (AGB) stars have been taken from van den Hoek & 
|Groenewegen| ( |1997| . 

Galaxies are identified as gravitationally bound groups 
of star particles by running a friends-of-friends (FOF) algo- 
rithm. Each FOF group is then decomposed into a set of dis- 
joint substructures, which are identified as locally overdense 
regions in the density field of the background main halo by 
the SUBFIND algorithm ( jSpringel et~aLl[200T] ) . After per- 
forming a gravitational unbinding procedure, only sub-halos 
with at least 20 bound particles are considered to be gen- 
uine structures ( |Saro et' al. 2006). For each galaxy in each 
of the snapshots of interest (z ss 5.7,6.6), we obtain the 
total halo/gas/stellar mass (M h /M g /M t ), the SFR (M»), 
the mass weighted gas/stellar metallicity (Z g /Z t ), the mass 
weighted age (t»), the mass weighted gas temperature and 
the half mass radius of the dark matter halo. 



2.1 Identifying LAEs 

The simulated properties of each galaxy at z as 5.7,6.6 are 
used to calculate the total intrinsic Lya (L™ ) and con- 
tinuum luminosity [L]. ) which include both the contri- 
bution from stellar sources and from the cooling of colli- 
sionally excited H 1 in the ISM as shown in |Dayal et al.| 
( |2010[ ). To calculate the stellar contribution, we obtain the 
spectrum of each LAE using the population synthesis code 
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STARBURST99 (Leitherer et al. 19991, using the simulated 



values of t„ , M* , M, ; the cooling radiation depends on the 
temperature of the ISM gas. The intrinsic Lya luminosity 
can be translated into the observed luminosity such that 
L a — L^faTa, while the observed continuum luminosity, 
L c is expressed as L c = L c nt f c . Here, f a (fc) are the frac- 
tions of Lya (continuum) photons escaping the galaxy, un- 
damped by the ISM dust and T a is the fraction of the Lya 
luminosity that is transmitted through the IGM, undamped 
by neutral hydrogen. 

The main features of the model used to calculate f c , f a 
and T a are: (a) for each galaxy the dust enrichment is de- 
rived by using its intrinsic properties (M„ t t , M g ) and as- 
suming Type II supernovae (SNII) to be the primary dust 
factories. The dust mass, Md, is calculated including dust 
production due to SNII, assuming an average dust mass pro- 
duced per SNII of 0.5 M dTodini & Ferrara||2001| iNozawa 



et al.||2003[ |2007| |Bianchi fc Schneider||2007^ dust destruc- 
tion with an efficiency of about 40% in the region shocked to 
speeds ^ 100 kms -1 by SNII shocks, assimilation of a ho- 
mogeneous mixture of dust and gas into subsequent SF (as- 
tration), and ejection of a homogeneous mixture of gas and 
dust from the galaxy due to SNII, (b) the dust distribution 
radius, r^, is taken to scale with the effective stellar distri- 
bution scale, r e , such that ~ (0.6, 1.0)r e at z ~ (5.7, 6.6) 
respectively; the calculation of r e is described later in Sec. 
3.1 (c) f c is calculated assuming a slab-like dust distribu- 
tion and we use f a = (1.5,0.6)/ c , as inferred for LAEs at 
z ~ (5.7,6.6), and (d) T a is calculated using the mean pho- 
toionization rate predicted by the Early Reionization Model 
(ERM, reionization ends at z ~ 7) as proposed by |Gallerani| 



et al. (20081, according to which the neutral hydrogen frac- 
tion X Hi = (6.0 x 10 _5 ,2.3 x 10" 4 ) at z « (5.7,6.6). 

Once the above calculations have been carried out, 
following the current observational criterion, galaxies with 
L a 10 42 ' 2 ergs _1 and an observed Lya equivalent width 
EW = L a /L c ^ 20 A are identified as LAEs. Complete de- 



tails of these calculations can be found in Dayal et al. ( 2010 1 



and Dayal & Ferrara (2011a I 



3 MOLECULAR HYDROGEN PHYSICS 

H2 can be formed in galaxies by two main methods: the first, 
and rather inefficient method involves gas-phase reactions 
mainly through the cou pled reactions e~ + H — > H~ + h v 
and H" + H -> H 2 + e" ( |McDowell|l96l"||Parla et al.|l983 ). 
The second, more efficient channel is through bond forma- 
tion on dust grains: this process begins with the collision 
and absorption of at least two hydrogen atoms by the same 
dust grain. The hydrogen atoms are weakly bound to the 
grain surface through Van der Waals forces and can mi- 
grate on the grain either by tunnelling or thermal hopping. 
If the hydrogen atoms encounter each other, bond formation 
takes place, the excess energy is released to the grain, and 
the H2 molecule is ejected into the gas phase ( Gould et al 
19631. However, H2 molecules so produced can be dissoci- 



ated by the far ultraviolet (FUV) interstellar radiation field 
in the Lyman- Werner (LW) band between 11.2-13.6 eV; the 
twin processes of self shielding and dust absorption (e.g. Hol- 
|lenbach fc Tielens|p99"l ) drive the shielding of H 2 to FUV 
photons, thereby preventing photodissociation. Considering 



these processes is of utmost importance since H2 is found in 
molecular clouds that are surrounded by a photodissociation 
region (PDR) where the gas is predominantly atomic. 

In this work we estimate the H2 mass of each LAE 
at z » (5.7, 6.6) using the analytic model presented in 
Krumholz et al.| ( [20081 [2009] ) and |McKee fc Krumholz| 
(20101, hereafter referred to as the KMT model. In brief, 



the KMT model considers an idealized spherical cloud im- 
mersed in a uniform, isotropic LW radiation field. Then, the 
equations of radiative transfer coupled to the H2 formation- 
dissociation balance are solved, assuming the cloud to be in 
the steady state. 

The analytical solution to the H2 mass fraction, /h 2 = 
Mh 2 /Mhi, is then obtained by solving for the radial position 
at which the transition between the atomic envelope and the 
molecular interior occurs within the cloud; in this equation 
Mh 2 and Mm refer to the mass of molecular hydrogen and 
the mass of neutral hydrogen respectively. The KMT study 
shows that the fraction of the radius at which this transi- 
tion occurs is solely a function of the dust optical depth in 
the LW band and the dimensionless parameter \, which are 
discussed in what follows. 



3.1 Modelling molecular hydrogen in LAEs 

We now describe the model used to calculate the H 2 con- 
tent of the LAEs identified in the simulation snapshots at 
z ~ 5.7, 6.6. We start by assuming that the MCs lie in a 
region that extends from the centre of the galaxy up to the 
effective stellar distribution radius, r e , calculated in |Dayai] 
(20101. This assumption has been motivated by the 



et al. 



fact that star formation occurs in MCs; the physical distri- 
bution scale of MCs and the stars is therefore expected to 
be quite similar. Further, the value of r e is based on esti- 
mates following the results of |Bolton et al.| (p008), who have 
derived fitting formulae relating the V-band luminosity and 
the stellar distribution scale from their observations of mas- 
sive, early type galaxies between z = 0.06 — 0.36. Though not 
an entirely robust estimate, we extend this result to galax- 
ies at z w 5.7 and 6.6 due to the paucity of observational 
data regarding the stellar distribution scales in high-redshift 
galaxies. However, such estimates are in surprisingly good 
agreement with recent observational results: Malhotra et al. 
( |2011[ ) find that the half-light radius of LAEs has a mean 
value ~ 0.16 arcsec at z ~ 5.7 and this remains constant 
for all redshifts in the range 2 ^ z ^ 6.5; these estimates lie 
within la of the mean value of the theoretical r e estimates 
used throughout this paper. 

As mentioned above, using the KMT model, the H2 
mass fraction solely depends on the dust optical depth in the 
LW band, t c , and the dimensionless parameter x, such that 
the analytical solution for the H2 fraction can be written as 
( |McKee fc Krumholz|20l"ol ): 



H 2 



1 



1 + 0.25s' 



(1) 



where the dimensionless parameter s can be expressed as 
_ ln(l + 0.6 X + 0.01x 2 ) 



0.6r c 



(2) 
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Table 1. As a function of the halo mass range (col. 1), we show 
the fraction of ISM gas with temperature T ^ 10 4 K (col. 2). 



M h 
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Wo] 




< 10 10 


0.58 


^qIO-10.4 


0.40 


IqIO.4-10.8 


0.35 


jqIO.8-11.2 


0.35 


> 10 112 


0.35 



Here, 



X 



fdissEpCTdC 



(3) 



where, fdiss —0.1 ( Draine fc Bertoldi|1996 1 is the fraction of 
absorbed LW band photons that lead to dissociation of the 
H2 molecules, Eq is the free space photon number density in 
the LW band, o~d is the dust absorption cross-section per hy- 
drogen nucleus to LW photons, c is the speed of light, «cnm 
is the number density of gas in the cold atomic medium that 
surrounds the molecular part of the cloud, and 1Z is the co- 
efficient for the rate of H2 formation on the surface of dust 
grains. 

It must be noted that these equations apply only for 
s < 2; for values of s ^ 2, the gas is predominantly atomic, 
such that /h 2 = 0. Also, we note that the calculations pre- 
sented here concern only average quantities in a spherically 
symmetric framework; a full calculation of the radial de- 
pendence of the model parameters is the subject of ongoing 
work. We now explain the calculations of Eo, o~d, bcnm, TZ 
and t c in what follows. 



3.1.1 LW photon number density 



As mentioned in Sec. |2.1[ we obtain the intrinsic spectrum of 
each LAE using the population synthesis code STARBURST99. 
Then, assuming all the stars to form at the centre of the 
galaxy, the number density of LW photons of a specific wave- 
length A (912 A < A < 1120 A), at a distance r from the 
centre can be expressed as 



n x (r) = 



4ircr 2 



(4) 



Here, L\ is the monochromatic luminosity at the wavelength 
A obtained using STARBURST99 and h is the Planck constant. 
The free space photon number density at radius r, in the 
entire LW band, ulw, can then be calculated by integrating 
over all the wavelengths in the band such that 



n L w{r) = 



1 

47rr 2 c 



1120A 



he 



dX. 



(5) 



The value of ulw averaged over a sphere of radius r e then 
gives the photon number density in the LW band such that 



E = (ntw) = 3n L w(r e ) 



(6) 



3.1.2 Cold neutral medium density 

Much of the neutral gas in galaxies is observed to be cold, 
with temperatures of order of ~ 100 K (CNM), or warm, 
with temperatures of order of ~ 10 4 K (WNM), in approx- 



imate pressure balance (Wolfire et al. 20031; as has been 
pointed out before, MCs form in regions where the gas is 
primarily cold. To calculate the density of the cold neutral 
medium, we start by obtaining the fraction, fj , of ISM gas 



with temperature T ^ 10 K (see also Fig. 1, Dayal et al 



20101; the value of f± averaged over galaxies of different 
halo masses is shown in Tab. [I] However, the large volume 
simulated (w 10 6 cMpc 3 ), naturally results in a low mass 
resolution, such that we are unable to resolve the cold and 
warm gas phases inside the ISM of individual galaxies. We 
therefore make the approximation that the ISM of each sim- 
ulated galaxy has an equal amount of cold and warm neutral 
gas. The mass of the cold neutral gas, Mcnm, in any galaxy 
can then be calculated as 



(7) 



where Mh = 0.76 M gas is the mass of hydrogen in the ISM. 
We assume the gas to be distributed in a disk with a radius 
r g and scale height H such that ( Ferrara et al.|2 000) 



H 



4.5Ar 2 oo, 
15.3A ( 

V Ue 



(8) 
(9) 



Here, value of the spin parameter is taken to be A = 0.04 



(Ferrara et al. 20001 and the virial radius, 7-200, is calcu- 
lated assuming the collapsed region has an overdensity of 
200 times the critical density at the redshift considered. As- 
suming that the typical velocity that determines the scale 
height of the disk is that of WNM, the effective gas sound 
speed is taken to be c s = lOkms -1 . Finally, the halo escape 
velocity, v e , is related to the circular velocity of the halo, v c> 
by the relation v e = 2pv c with p — 1.65 (Mac Low & Fer- 
rara 1999). The average global number density of the CNM, 
ncNM can then be expressed as 



icnm — 



Mqnm 
■nriHmH ' 



(10) 



where ran is the hydrogen mass. 



3.1.3 Molecular hydrogen formation rate 

The rate of H2 formation on dust grains, 1Z, can be expressed 



as ( Hirashita & Ferrara||2005 1 



n 



■- 4.1 x 10~ L 'S 
V100K 



0.1 /im 

S 

2 gem -3 



V 

Trp 



(ii) 



where 5* = S(T, Td) is the sticking coefficient for hydrogen 
atoms on dust grains, T is the gas temperature, 8 is the 
density of the dust grains, T> is the dust-to-gas ratio, and 
a is the radius of the dust grain which is assumed to be 
a sphere. Following the assumptions of |Dayal et al.| ( [2010[ ), 
we assume all the dust grains to be carbonaceous such that 
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5 — 2 gem 3 ; using the size distribution of SNII-produced 
dust grains ( |Todini fc Ferrara|2001[ ), we use an average gram 
size value of a = 300 A. As mentioned before, dust predom- 
inantly forms in high-density, cold MCs which have a more 
effective self-shielding to LW photons compared to the lower 
density WMN; following this argument, we use T = 100 K. 



Further, the sticking coefficient S(T, T d ) is given by Hoi 



lenbach & McKee ( 1979 \ and Omukai ( 2000 1 as 



S(T, T d ) = [1 + 0.04(T + T d ) ' 5 + 2 x 10~ 3 T + 8 x 10~ 6 T 2 ] -1 

x{l + ea;p[7.5 x 10 2 (l/75 - l/T d )]}-\ 

(12) 

where T d is the dust temperature. However, the sticking co- 
efficient is not affected by the exact dust temperature as 
long as T d sC 70 K, which is true for all LAEs, which have 
dust temperatures T d w 30 — 45 K. 

3.1.4 Dust cross-section and optical depth 

We now discuss the calculation of the dust cross-section per 
hydrogen atom to LW photons, which can be expressed as 



a d = 



(13) 



where a = 300 A is the average radius of SNH produced dust 
grains as mentioned above, Q a bs{a) is the effective cross- 
section normalised to the geometric one for the average grain 
size and a wavelength of 1000 A corresponding to the centre 
of the LW band ( Draine fc Lee|1984 | and Nh is the number 
of hydrogen nuclei per dust grain. To calculate Nh, we use 
the dust to gas ratio, T> such that 



V : 



Ua 3 SN dus 



M dua 



M gaa fim p N H 



(14) 



where /1 = 0.59 is the mean molecular weight of a fully 
ionized gas of primordial composition. Substituting Nh from 
Eq. [13] into Eq. [H] yields 

= — ■ (15) 

From Eqns. [IT] and [15) we note that x i s independent of T>. 

Finally, the dust optical depth in the LW band is cal- 
culated as: 

_ J2 d Qabs(a) 



i a5 



(16) 



where the dust surface density T, d = M d /(irr d ) and r d « 
(0.6,1.0) x r e at z w (5.7,6.6) as required to best fit the 
LAE U V luminosity functions ( Dayal et al.|2010 1 . 

Once these calculations have been carried out, the total 
molecular hydrogen mass, Mh 2 in each LAE is estimated as 



Mh 2 = /h 2 M H i, 



(17) 



where Mm = f^Mn is the neutral hydrogen mass in the 
gas disk within a radius r e and scale height H, where star 
formation takes place. 



4 RESULTS 

Once that the above calculations have been carried out, wc 
can discuss the results regarding the molecular hydrogen 



fraction and the total H2 mass, and relate these to the phys- 
ical properties of LAEs at z ~ 5.7,6.6. In what follows, we 
also calculate the visibility of such LAEs through their CO 
emission and end by making predictions for such detections 
in the ALMA Early Science Release (ESR) . 

4.1 Molecular hydrogen content of LAEs 

As shown in Sec. |3.1| the H2 fraction is decided by three 
important physical parameters: (a) the SFR which deter- 
mines the intensity of the H2-dissociating LW field, (b) the 
amount of cold gas available to shield the MC against the 
LW field, and (c) the dust mass on which the H2 forms, 
and which additionally shields the molecular hydrogen by 
absorbing LW photons. We now quantify how the molecular 
hydrogen fraction depends on each of these parameters. 

We start our discussion by mentioning that the stel- 
lar mass of LAEs ranges between M» = 10 8 ~ 10 - 4 M Q at 
z ~ 5.7. As expected in a hierarchical structure forma- 
tion scenario where progressively larger objects form from 
the merger of smaller ones, such range narrows to M* « 
1O 81 " 1O M at z w 6.6 (see also Tab.Sh; galaxies with stel- 
lar masses above 1O 1O M0 have not had the time to assemble 
in large numbers by z w 6.6. Further, the SFR of LAEs falls 
in the interval M* » 0.8 - 12OM yr~ 1 at z » 5.7, with 
larger (M* ^ 1O 9 ' 5 M0) galaxies being the most efficient in 
star formation. At stellar masses lower than this value, the 
SFR-stellar mass relation flattens at both the redshifts con- 
sidered, as seen from Fig. [l] due to the stronger effects of 
mechanical feedback (ejection of gas in outflows) inhibiting 
star formation (see |Dayal et al^]p00 9). As a result of the 
narrower stellar mass range, the SFR for z m 6.6 are con- 
centrated in a narrower range between 1.6 — 46 Mq yr~ (see 
Tab. [2}. 

As for the molecular hydrogen fraction, starting with 
z « 5.7, we find /h 2 ^ 0.1 for galaxies with M, < 10 9 M Q ; 
it rises to ~ 0.5 — 0.9 for intermediate mass galaxies with 
M* = 10 9_10 M Q and then decreases again to 0.2-0.6 for 
the few largest galaxies. Such behavior can be explained 
as follows: although the smallest galaxies (M» < 10 9 Mq) 
have the smallest SFR (and hence weakest LW field), they 
also less dusty, resulting in a lower H2 production rate and 
self-shielding ability against photodissociation. On the other 
hand, the larger SFR compared to the dust and cold gas 
mass in the largest galaxies leads to a decreased /h 2 ■ It is 
thus the intermediate mass galaxies that show the largest 
H2 fraction by virtue of a well-tuned balance between the 
H2 formation and dissociation rates. Such argument is sup- 
ported by the fact that at a given value of the stellar mass, 
intermediate mass galaxies with the lowest SFR have the 
largest H2 fraction as seen from Fig. [l] Qualitatively, the 
situation remains the same at z w 6.6, although quantita- 
tively, the /h 2 value is lower for all LAEs; the reason for this 
is detailed in what follows. 

To understand the /h 2 dependence on the cold gas mass 
(Mcnm), we start by noting that Mqnm scales well with the 
SFR as shown in Fig. [2] more massive galaxies have a smaller 
cold mass fraction, possibly due to stellar sources heating a 
greater part of the ISM to higher temperatures (see also Tab. 
[T]). From the same figure, we see that for a given value of 
Mcnm, galaxies with the lowest SFR have the largest value 
of /h 2 which is easy to understand considering that for a 
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Table 2. For all the LAEs at the redshifts shown (col. 1), we show the range of stellar mass (col. 2), the range of SFR (col. 3), the 
range of cold neutral medium mass (col. 4), the range of dust mass (col. 5), the average molecular fraction (col. 6), the average mass of 
molecular hydrogen (col. 7), the average intrinsic CO(l-O) luminosity (col. 8), and the average value of CO(6-5) flux (col. 9). 
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Figure 1. The molecular hydrogen fraction, /h 2 i of LAEs at 
z ft: 5.7 (6.6) as a function of the stellar mass (M„) and SFR 
(M„) are shown in the upper (lower) panels respectively. Points 
are color-coded for different values of /jj., • 



Figure 2. The molecular hydrogen fraction, /h 2 , of LAEs at 
z RJ 5.7 (6.6) as a function of the stellar mass (M*) and mass 
of cold H I ISM gas with temperature less than 100 K (Mcnm) 
are shown in the upper (lower) panels respectively. Points are 
color-coded for different values of /jj 2 ■ 



given Mcnm, at lower SFR, the H2 LW dissociation becomes 
less efficient. Again, it is the intermediate mass galaxies that 
have the largest value of /h 2 , while the H2 fraction is larger 
at z « 5.7 compared to z « 6.6 for all LAEs. 

Dust also plays a key role in terms of the H2 abundance 
since this molecule predominantly forms on dust grains 
which also shield the molecular hydrogen so formed, by ab- 
sorbing LW photons. The dust mass for each LAE has been 
calculated as explained in Sec. |2.1[ since we assume stellar 
sources (i.e. SNII) to be the main dust producers, the total 
dust mass scales well with the stellar mass, ranging between 
10 3 ' 4_7 ' 2 Mq at z « 5.7. As mentioned before, since H2 forms 
on dust grains, naively it might be expected that the larger 
the dust content, the larger the value of /h 2 - However, this 
does not hold true; the LW intensity in the largest galaxies 



is enough to dissociate the H2 formed as a result of which, 
the intermediate mass galaxies end up with the largest H2 
fraction. As expected, at a given value of M*, galaxies with 
the largest dust mass have the largest /h 2 value (Fig. |3p . 

Finally, we discuss why the value of /h 2 is lower at 
z w 6.6 compared to z r; 5.7, even though galaxies iden- 
tified as LAEs are extremely similar at these redshifts, in 
terms of M», M», Md and A/cnm (see Tab. [2I: as men- 
tioned in Sec. |3.1.4[ the value of the dust distribution ra- 
dius, rd = (0.6, 1.0)r e at z m (5.7, 6.6), is fixed by matching 
the observed LAE UV luminosity functions. These values 
imply that the dust is more concentrated in the inner parts 
of LAEs than the stars themselves, perhaps hinting at the 
existence of dust/ met allicity radial gradients. We note that 



Molecular hydrogen in LAEs 7 



/h 2 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 




1.5 9.0 9.5 

l_og(M„ / M a ) 



10.0 



10.5 



Figure 3. The molecular hydrogen fraction, /h 2 i of LAEs at 
z fa 5.7 (6.6) as a function of the stellar mass (M t ) and dust 
mass (M(j) are shown in the upper (lower) panels respectively. 
Points are color-coded for different values of fn . 



r e , the stellar distribution scale, is similar at both the red- 
shifts considered. The larger value of r& results in a dust 
optical depth that is about 1.6 times smaller at z ~ 6.6 as 



compared to that 5.7 (see Eq. 161; although comparable H 



masses would be produced in LAEs with similar physical 
properties at both z fa 5.7 and 6.6, a larger amount is disso- 
ciated at z ~ 6.6 due to a decreased dust absorption of HV 
dissociating LW photons. Averaged over all LAEs, /h 2 ~ 0.6 
at z fa 5.7, and only /h 2 ~ 0.3 at z fa 6.6, as seen from Fig. 
[4] and Tab. [2] From the same figure, we see that while /h 2 
covers the broad range — 0.85 at z fa 5.7, no LAEs have 
/h 2 0.65 at z fa 6.6, as a result of the smaller dust optical 
dept h. F inally, translating the value of fu 2 into a total Mh 2 
(Eq 




Figure 4. Normalized distribution of the number of LAEs as 
a function of the molecular hydrogen fraction /n 2 at z ~ 5.7 
(dark/blue) and z fa 6.6 (light/red). 



4.2 CO detectability in LAEs 

Now that the dependence of /h 2 on the physical properties 
of LAEs has been understood, we are in a position to make 
predictions for the H2 detectability in these galaxies. H2 
has strongly forbidden rotational transitions; the rotational- 
vibrational lines have very high excitation temperatures, 
that are attainable only under somewhat extreme conditions 



involving intense irradiation or shock waves (see Solomon <fc 
|Vanden Bout|2005[ ) in the absence of which the H2 is invis- 
ible. On the other hand, CO has a weak dipole moment: its 
rotational levels are then easily excited and thermalized by 
collisions with H2. In addition, CO is a very stable molecule 
and the most abundant after H2 (e.g. Solomon & Vanden 



Bout|2 005 Om ont|2007 |. Because of such considerations, it 
is popularly used as a tracer of H2. 

The luminosity of the CO(l-O) transition can be related 
to Mh, mass as 



Leo = Mn 2 /a. 



(18) 



The parameter a used in the equation above depends on the 
distribution of star forming clouds: in the Milky Way, where 
star formation takes place in molecular clouds with dense 
cores, confined by self gravity, 01 = 4.6 Mq K kms -1 pc 2 
| |Solomon et al.|[l987[ ). On the other hand, in high-redshift 
Ultra Luminous Infrared Galaxies (ULIRGs), where star for- 
mation is expected to occur in a dense intercloud medium 
bound by the potential of th e galaxy, a has a much lo wer 



value of 0.8M s Kkms" 1 pc i (Downes & Solomon 1998 1. It 



5.7 and 6.6 respectively, as shown in Tab. [2] 
To summarize, we find that intermediate mass LAEs 
(M, fa 1O 9_1O M0) have the largest molecular hydrogen frac- 
tion; a delicate balance between M» , Mcnm and ALj leads to 
H2 formation (and self-shielding) dominating over H2 photo- 
dissociation. Finally, due to a smaller dust optical depth (by 
a factor fa 1.6) at z fa 6.6, the H2 fraction is about twice as 
large at z fa 5.7 for LAEs with comparable physical proper- 
ties. 



is worth noting that the value of a is also related to the 
metallicity of the inte rstellar gas (e.g. |Leroy et al.||2009| 



17 1, we find the average value of Mh 2 fa 10 8 ' 9 , 10 8 ' 4 at Narayanan et al.|2011 Genzel et al.|2011 1. Using the Green 



Bank telescope, Wagg et al. (2009) undertook a search for 



CO emission in two LAEs at z > 6, and they adopted the 
ULIRG value of the conversion factor to estimate Mh 2 . For a 
reasonable comparison with such LAE data, we use the same 
value of a = 0.8 M© Kkms -1 pc 2 to compute the CO(l- 
0) luminosity values for simulated LAEs at z fa (5.7, 6.6). 
We find that the value of Leo scales with L a for both the 
redshifts considered as shown in Fig. [5] quantitatively we 
find, Leo oc La' 04 ' 1 at z fa (5.7, 6.6) respectively. Al- 
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Figure 5. As a function of the observed Lyct luminosity (L a ), we 
show the estimated CO luminosity at z fti 5.7 and 6.6 in the upper 
and lower panels respectively; the points shown are color-coded 
according to the SFR. In the lower panel, the arrows indicate 
the upper limit on the L C o obtained by |Wagg et al.| | |2009| l: (i) 
corresponds to the nondetection of the CO ( 1-0) line of HCM 6A, 
a LAE at z ft 6.6, with SFR ft 2 M0yr _1 | |Hu et al.|2002| . (ii) 
corresponds to the nondetection of t he CO (1-0) line of IOK-1, a 
z ft 7 LAE with SFR « 10 M yr- 1 jlye et al.|2006| ). 



though such a relation has a huge scatter, this implies that 
the brightest LAEs are the best candidates for molecular 
emission searches. 

Interestingly, such a trend also implies that although 
the value of /h 2 is the largest for the intermediate 



mass/luminosity LAEs (see Sec. 4.1|, the increasingly large 
neutral hydrogen mass, A/hi, in the star forming disk (see 
Eq. 17 1 of massive galaxies wipes out such a subtle signa- 
ture, whose imprint remains as a flattening of the Leo — L a 
relation towards the most luminous objects. As expected, as 
a result of their lower dust optical depth, and hence, a lower 
molecular hydrogen fraction, the CO luminosity for LAEs is 
smaller at z ~ 6.6 than at z ~ 5.7; averaged over all LAEs, 
the CO luminosity at z ~ 6.6 is about a factor 3 lower than 
that at z « 5.7, as shown in Tab. [2] As a validation of our 
model, our theoretical Leo estimates are in accord with the 
upper limits found by Wagg et al.| ( |2009 \ for their observed 
LAEs, as seen from lower panel of Fig. |5j which represents 
an encouraging sanity check of our model. Finally, we note 
that a is expected to be larger than the Galactic value in 
low metallicity environments ( Leroy et al.|20 09 Narayanan 



|et al.|[20TT] |Genzel et al~1|2011[ ), which is likely the case for 
LAEs, while we have used a value derived using ULIRGs 



As expected, an increase in the value of a would lead to a 
decrease in the CO luminosity, thereby negatively affecting 
the CO detectability of LAEs. 

4.3 Predictions for ALMA early science 

Now that the CO(1-0) luminosity has been calculated for 
all LAEs in our simulation, we make predictions for the de- 
tectability of such CO lines using ALMA. As of now, these 
observations have been limited to the most luminous high-z 



sources such as QSOs (e.g. 


Cox et al.||2002 


Bertoldi et al. 


2003; Walter et al.|2004; Wcifi ct al. 2007||Wang et al. 2010 


Riechers et al. |2011|) and sub-millimeter g 


alaxies (SMG) 
normal galax- 


(e.g. Greve et al.||2005 Tacconi et al.||2006l; 



ies such as Lyman break galaxies (LBGs) have only been 
detected at a much lower redshifts, z w 3. We start by 
noting that ALMA Early Science (Cycle 0) consists of the 
use of 16 antennas and a limited number of Bands. In this 
configuration, the lowest CO rotational transition observ- 
able at z « 5.7 is CO(5-4) {u nst = 576.267GHz), which 
falls within ALMA Band 3 (84 - 116 GHz). However, as is 
shown in what follows, the rotational transition strength of 
the CO (5-4) line is quite comparable to that of CO (6-5). 
The latter is the lowest line observable with ALMA Cycle 
for both (z ~ 5.7, 6.6) samples, for which reason we show 
results for the CO(6-5) transition in what follows; in the 
future, using the the full capabilities of ALMA, lower fre- 
quency bands (31.3 — 45, GHz and 67 — 90 GHz) will also 
allow for the study of lower-J CO line transitions at high 
redshift. 

We now describe how we translate the CO(1-0) luminosi- 



ties calculated above in Sec. 4.2 in to CO(6-5) luminosities 



Theoretical fits to observational CO spectral energy distri- 
butions (SEDs) have been carried out by several authors 
using Local Thermal Equilibrium (LTE; Obreschkow et al. 
2009 p or Large Velocity Gradient (LVG) models ( Bayet et al. 



2009 p to describe the molecular gas. We use the model pro- 
posed by Obreschkow et al. (2009) which assumes a single 



gas component in LTE. In this framework, the frequency in- 
tegrated luminosity from the transition ( J — > J — 1) can be 
expressed as ( Obreschkow et al.|2009 l: 



L( J — > J — 1) oc 1 — exp(rj) 



J 



(19) 



exp 



where T e — 100 K is the gas excitation temperature, z/co = 
115.271 GHz is the rest-frame frequency of the CO(1-0) tran- 
sition, t.j is the optical depth and k b is the Boltzmann con- 
stant. Further, tj can be expressed as 



t,j = 7.2rexp(- 



hvcoJ^ 
2k b T e 



)sinh 



hvcoJ 
2k b T e 



(20) 



where r is an experimental determined normalization con- 
stant, which we take to be f = 2, following the results of 
breschkow et al.| ( (2009| ). 

Considering that the frequency of the ( J — > J — 1) tran- 
sition is related to that of (1-0) by vj — JVco, and using 
the relation between frequency-integrated luminosity L, and 
the brightness temperature luminosity Leo (cfr. App. A, 
Obreschkow et al.|2009| ): 



Leo = (87rftb) A e L = (S-nkbY 



(21) 
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Figure 6. The line integrated CO(6-5) flux, Sco, as a function of 
the L a observed for LAEs at z ~ 5.7 (circles) and z 6.6 (trian- 
gles). Horizontal lines represent the ALMA Cycle sensitivity for 
a signal-to-noise ratio of 5, for an integration time of 1 hour (dot- 
dashed line), 5 hours (dotted line) and 10 hours (dashed line), 
assuming a spectral resolution of 50 km/s. Solid line represents 
the sensitivity (S/N = 5) achieved with 10 hours of integration 
time with 50 antennas. 



we obtain the Leo (6-5) luminosity in units of [K kms 1 pc 2 ] 
and find that L CO (6-5)= 1.58 Leo (1-0), while L CO (5-4)= 
1.75 Leo (1-0); all the results presented for the C0(6-5) tran- 
sition also remain largely valid for the CO(5-4) transition. 

This line luminosity can be converted into the line in- 
tegrated flux, Sco such that 



L co = 3.25 x 10 7 ScoAV*/- 6 2 s (l + z)~ 3 D 2 L 



(22) 



where Sco&V = S v is the velocity integrated flux 
[Jykms -1 ], v bs is the observed frequency in GHz, and Dl 
is the luminosity distance. We then assume that the CO(6-5) 
line has a gaussian profile with a width given by the rota- 
tional velocity of the galaxy, v r , which we take to be equal to 
1.5 times the halo rotation velocity (see Dayal et al.|2009 |. 
The average value of v r for LAEs at both the redshifts con- 
sidered is of order of 200 km/s. The observable CO flux Sco 
can then be expressed as Sco — S v /v r . 

As noted in Sec. |4.2| the CO(1-0) luminosity value scales 
with the observed Lya luminosity; the line integrated flux 
is therefore also expected to behave in a similar way. This is 
indeed the case, as seen from Fig. [6] Albeit with a large scat- 
ter, galaxies with the largest Lya luminosity show the largest 
value of Sco (6 — 5) at both z w 5.7, 6.6; as expected from a 
comparison of the CO(1-0) luminosities, the average value of 
Sco is about 3 times lower at z ~ 6.6, compared to z w 5.7 
(Tab.[2|. We find that atz« 6.6 none of the LAEs would be 
detectable in CO, even with an ALMA Cycle integration 
time for about 10 hours. On the other hand, at z ~ 5.7, 
about 1-2% of LAEs, i.e., those with L a > 10 43 2 ergs _1 
could be detectable with an ALMA integration time of 5-10 
hours (i.e. a detection limit of ~ 1-1.4 mjy respectively), 
assuming a signal-to-noise ratio S/N = 5. However, using 
50 antennas of ALMA with an integration time of 10 hours 
increases the sensitivity to about 0.3 mjy making a signif- 



Figure 7. The number density of LAEs at z ss 5.7 and 6.6 are 
shown as a function of the line integrated CO (6-5) flux using 
dashed and solid lines respectively; shaded regions show the pois- 
sonian errors. Vertical dot-dashed, dotted and dashed lines rep- 
resent the ALMA sensitivity limits for 16 antennas assuming a 
signal to noise ratio, S/N = 5, for an integration time of 1, 5 and 
10 hours respectively. Vertical solid line represents the sensitivity 
(S/N = 5) achieved with 10 hours of integration time with 50 
antennas. 



icant change such that about 13%, 1.4% of LAEs become 
detectable at z « 5.7 and 6.6 respectively. 

To clarify such results, we present the CO 'flux func- 
tion', i.e. the number density of LAEs as a function of the 
line integrated CO flux, Sco (6 — 5). As is clearly seen again, 
none of the z fa 6.6 LAEs are detectable with the Cycle 
of ALMA, even for a 10 hour detection limit of about 1 
mjy. On the other hand, the z « 5.7 flux function, extends 
into the 1-10 ALMA detection bands as mentioned above. 
However, the peak of the function lies at Sco ~ 0.1 mjy at 
z ss 5.7; due to the lower H2 mass at 2 « 6.6, the flux func- 
tion peaks at a lower value of about 0.06 mjy. As mentioned 
above, using 50 ALMA antennas for an integration time of 
10 hours increases the number of LAEs detectable in CO; 
as expected, this leads to a larger part of the 'flux function' 
being observed at both the redshifts considered. 



5 SUMMARY AND CONCLUSIONS 

Using a physically motivated model, we present results con- 
cerning the H2 fraction, and the total H2 mass inside the 
ISM of LAEs at z w 5.7 and 6.6 by coupling a semi-analytic 
model of the molecular content of galaxies, to a previously 
developed LAE model that reproduces a number of observed 
LAE data sets. 

We start by using cosmological SPH simulations to ob- 
tain the physical properties of each simulated galaxy, includ- 
ing the total halo/gas/stellar mass (M h /M g /M*), the SFR 
(M»), the mass weighted gas/stellar metallicity (Z g /Z t ), the 
mass weighted age (t t ) and the mass weighted gas tempera- 
ture; the values of Z m ,t», M, of each galaxy are then used to 
produce its intrinsic spectrum using STARBURST99 ( Leithcrcr 
et al. 19991. Assuming SNII to be the primary sources of 
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dust, we calculate the ISM dust attenuation of both Lya 
and continuum photons, and the IGM Lya transmission, to 



select galaxies as LAEs (Dayal et al. 2009 2010 Dayal & 



Ferrara 2011b I. For each of these galaxies, the value of /h 2 is 



calculated using the analytic model proposed by |Krumholz| 
et al.| (|2008[ |2009|) and |McKee fc Krumholz| (|2010|, which 



considers an idealized spherical molecular cloud immersed 
in a uniform and isotropic H2-dissociating Lyman- Werner 
radiation field; under approximations of steady state, the 
value of /h 2 is solely a function of the dust optical depth to 
LW photons, and a dimensionless parameter, \- The latter 
depends on: (a) the intensity of the LW field which is deter- 
mined by AL», (b) the total dust mass, Md, which both en- 
ables H2 formation on the surface of dust grains, and shields 
the H2 to LW photons, and (c) the cold H 1 mass, Mcnm, 
providing the key ingredient to produce H2. 

Using this model, we find that at z ~ 5.7, the value 
of /h 2 peaks and ranges between 0.5 — 0.9 for intermedi- 



ate mass LAEs with M, 



10 M 



Mq, decreasing for both 



smaller and larger galaxies; this trend also holds at z ~ 6.6. 
Such behaviour can be explained as follows: compared to in- 
termediate mass galaxies, smaller galaxies have smaller SFR 
(and hence a lower value of the dissociating LW field), lower 
dust masses for forming H2, and a lower cold gas (and dust) 
mass for self-shielding the H2 so formed against the dissoci- 
ating LW field; on the other hand, the larger SFR compared 
to the dust and cold gas mass in the largest galaxies leads 
to a lower /h 2 - Such an argument is validated by the fact 
that at a given stellar mass, galaxies with the lowest SFR, 
largest cold gas mass, and largest dust mass have the largest 
value of /h 2 . Further, we find that for LAEs with compara- 
ble SFR/gas mass/dust masses, /h 2 is about twice as large 
at z m 5.7 than at 6.6; to reproduce the LAE UV luminos- 
ity function data, the dust distribution radius at z « 6.6 
is about 1.6 times larger than that at 5.7, for similar LAE 
properties. As a result, the dust optical depth to LW pho- 
tons is lower at z ss 6.6 which leads to a larger amount of 
H2 being dissociated. 

We then translate the H2 mass we obtain into a CO 



luminosity. Recently, Wagg et al. (20091 used the green 



bank telescope to look for CO emission in two LAEs at 
z > 6 both of which resulted in non-detections. To com- 
pare our model predictions to their observations, we used a 
value of a = 0.8 Mq Kkms -1 pc 2 to translate the H2 mass 
into a CO(1-0) luminosity, Leo- We find that the value of 
Leo scales with L a , for both the redshifts considered; al- 
though the value of fui is the largest for the intermediate 



mass/luminosity LAEs (see Sec. 4.1|, the increasingly large 



neutral hydrogen mass, Mhi, in the star forming disk of 
increasingly massive galaxies wipes out such a subtle signa- 
ture. This result also implies that the brightest LAEs are 
the best candidates for molecular emission searches. As a 
reasonable validation of our model, our theoretical Leo es- 
timates are in accord with the upper limits found by | Wagg] 
et al. ( 2009 I in their experimental work. 

At z ~ 5.7, 6.6, the lowest CO rotational transition ob- 
servable with ALMA is the CO(6-5). We find that at z ^ 5.7, 
about 1-2% of the LAEs, i.e. those with L a > 10 43 2 ergs _1 , 
could be detectable with an integration time of 5-10 hours 
respectively, assuming a signal-to-noise ratio S/N = 5. Our 
results at z ~ 6.6 are more pessimistic; none of the LAEs 
would be detectable in CO, even with an ALMA integra- 



tion time of about 10 hours. We also present the CO 'flux 
function', the number density of LAEs as a function of 5co 
where we show that the number density of objects peaks at 
a value of about 0.1 mjy, which is much beyond the sensitiv- 
ity of ALMA; this peak shifts to progressively lower values 
with increasing redshift. 

Finally, we discuss the main caveats in the model. First, 
the calculations presented here concern only average quanti- 
ties in a spherically symmetric framework and a full calcula- 
tion of the radial dependence of /h 2 is the subject of an ongo- 
ing work. Secondly, the dust masses used in this work have 
been calculated assuming SNII to be the primary sources 
of dust production which is a reasonable assumption given 



that a number of authors (Todini & Ferrara 2001 Dwek 



et al.[2 007) have shown that the contribution of AGB stars 



becomes progressively less important and at some point neg- 
ligible towards higher redshifts (z ^ 5.7) when the Universe 
is less than 1 Gyr old. However, it must be noted that under 
certain conditions thought to hold in quasars, in which ex- 
tremely massive starbursts occur, the contribution of AGB 
can become important somewhat earlier (see Valiante et al. 
|2009[ > . Thirdly, the stellar distribution scale (and radial ex- 
tent of the MC) is based on estimates following the results 
of Bolton et al. (20081, who have derived fitting formulae 



relating the V-band luminosity and the stellar distribution 
scale from their observations of massive, early type galaxies 
between z — 0.06 — 0.36; however, we note that such esti- 
mates are in surprisingly good agreement (to within I— a) 
with recent observational results of z w 5.7 LAEs (Malhotra 
et al.|[2011| . Fourthly, the large cosmological volume simu- 



lated naturally results in a low mass resolution, such that 
we are unable to resolve density and temperature of the gas 
in the interior of individual galaxies. We have therefore esti- 
mated the ISM gas temperature distribution over broad halo 
mass bins. Further, in our calculations, we have assumed 
that half of the gas mass with temperatures of T < 10 4 K is 
cold, with T < 100 K. Finally, we have used a ULIRG value 
for the factor a, used to translate the H2 mass into a CO 
luminosity. However, this parameter depends on a number 
of poorly known properties of high- z galaxies, such as the 
spatial/mass distribution of their MCs and the gas metallic- 
ity ( Leroy et al.||200"9 Narayanan et al7||201 1 Genzel et al 



2011 1. It is hoped that upcoming data from state of the art 



instruments such as ALMA will be able to shed light and 
clarify such thorny issues. 
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